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Abstract 

We construct explicit Darboux transformations for a generalized, two-dimensional 
Dirac equation. Our results contain former findings for the one-dimensional, stationary 
Dirac equation, as well as for the fully time-dependent case in (1+1) dimensions. We 
show that our Darboux transformations are applicable to the two-dimensional Dirac 
equation in cylindrical coordinates and give several examples. 



1 Introduction 

The famous Darboux transformation [1] [2] is a mathematical scheme that interrelates 
solutions of differential equations by means of a linear differential operator. Since the 
Darboux transformation does not involve changes of variables, it is essentially different 
from other solution-generating methods, such as approaches related to Lie symmetries. 
While the original Darboux transformation applied to linear, one-dimensional equations 
of second order only, in the meantime it has been generalized to a variety of linear and 
nonlinear models in higher dimensions [4] [5]. While such nonlinear equations include 
the nonlinear Schrodinger equation, the nonlocal Korteweg-de-Vries equation or the sine- 
Gordon equation, in the context of linear models Schrodinger-type equations and the Dirac 
equation are particularly important. The Darboux transformation has been extensively 
applied to Schrodinger equations, especially since it has been found to be equivalent to the 
quantum-mechanical supersymmety formalism, see [3] and references therein. As far as the 
Dirac equation is concerned, a Darboux transformation has been constructed for both the 
one-dimensional stationary case [7] and the fully time-dependent equation [6] [9] in (1+1) 
dimensions. A further generalization of the Darboux transformation to Dirac equations in 
higher spatial dimensions will be made in the present work. In particular, we focus on a 
generalized Dirac equation in two variables, a special case of which is the stationary, two- 
dimensional Dirac equation. The Darboux transformation for our generalized equation 
will be constructed by means of an approproate intertwining relation (section 2), such 
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that the transformed solution as well as the transformed potential emerge in explicit form. 
A particularly interesting application of our Darboux transformation concerns the two- 
dimensional Dirac equation in cylindrical coordinates. This equation will be derived in 
section 3 and shown to admit our Darboux transformation in each of its variables (section 
4). In section 5 we state two examples of how the Darboux transformation can be applied 
to the two-dimensional Dirac equation in cylindrical coordinates. Starting from the force- 
free case, we obtain a Dirac equation minimally coupled to a magnetic field, together with 
a position-dependent mass. 

2 Construction of the Darboux transformation 

In the following we will introduce our generalized Dirac equation and construct a Darboux 
transformation by means of the intertwining relation. 

The Dirac equation. The two-dimensional equation we will study in this work has the 
following form: 

(M d x + N d y + Wo) V = 0, (1) 

where the symbol d denotes partial differentiation. The function i/j = ifj(x,y) stands for 
the 4-component solution vector and M = M(x,y), N = N(x,y), and Wo = Wo(x,y) 
denote 4x4 matrix functions. Furthermore, the function Wq will have the usual form 
Wo = E — Vo, where E is the constant energy and Vo = Vo(x, y) represents the potential. 
Observe that the (1+1) dimensional Dirac equation is a special case of (P), if we set 

E = M = - iffl (° *) (2, 

where a%, o<i stand for the Pauli matrices, and I2 denotes the 2x2 identity matrix. 
Furthermore, under this reduction y assumes the role of the time variable. 

The Darboux transformation. We will now construct the Darboux transformation 
for our Dirac equation ([1]). To this end, let us first multiply the equation by iV -1 on both 
sides, which will simplify our subsequent calculations. After the multiplication, equation 
(HJ takes the form 

(F 8 x + d y + Uo) v = 0, (3) 

where F = N~ 1 M and Uo = N^ 1 Wo- Now consider another Dirac equation of the form 
d3]), but for a different potential term U\ = U\{x,y): 

(f d x + d y + t/i) = 0, (4) 

where <f> = <j)(x, y) stands for the solution. Our goal is to map solutions of Q onto solutions 
of dH) by means of a linear, first-order differential operator L: 

L = A d x + A B. (5) 

Here the functions A = A(x,y) and B = B(x,y) must be determined in such a way that 
the solutions ip and 4> of the Dirac equations ([3]) and respectively, are related to each 
other via 

cp = Lip. 
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This is true if L satisfies the following operator equation, called intertwining relation: 

(Fd x + d y + U 1 )L = L (F 8 X + d y + Uq) . (6) 

A solution L of this relation is called a Darboux operator or intertwiner for equations (|3j) 
and ([3]). Now, in order to find an intertwiner of the form ([5]), we substitute its explicit 
form into relation ([6]) and expand both sides. Starting with the left hand side, we obtain 
the following expression: 

(F d x + 3 y + U x ) L = (Fd x + d y + Ux) (A d x + A B) 

= F A d xx + A d xy + {F A x + F A B + A y + Ux A) d x + 
+ A B d y + F A x B + F A B x + A y B + A B y + Ux A B. (7) 

In the same fashion we expand the right hand side of ©: 

L (F d x + d y + Uq) = {A d x + A B) (F 8 X + d y + Uq) 

= A F d xx + A d xy + (A F x + A Uq + A B F) d x + 

+ ABdy + A (U ) x + ABU . (8) 

Our intertwining relation ([6]) is fulfilled if the expressions on its left hand side ([7]) and on 
its right hand side (|8|) are the same, that is, if the coefficients of the respective derivative 
operators are equal on both sides. Comparison of these coefficients gives the following 



conditions: 

FA = A F (9) 

F A x + F A B + Ay + Ux A = AF X + AU + ABF (10) 

F A x B + F AB X + A y B + ABy + UxAB = A(U ) X + ABU . (11) 



We will now solve these equations for the functions A, B and for the transformed potential 
Ux- The first equation ([9]) is fulfilled if A and F commute, a fact that we will use for the 
simplification of terms in subsequent calculations. The second equation ([10]) can be solved 
with respect to the transformed potential Ux- 

Ux = A(F X + U + B F - F A' 1 A x - F B - A' 1 Ay) A' 1 . (12) 

It remains to solve the last condition (|lip . To this end, we first redefine the function B 
that appears in the intertwiner ([5]) as follows: 

B = -u x u-\ (13) 

introducing a 4 x 4 matrix u = u(x,y), the meaning of which will be determined below. 
Next, we insert the latter redefinition of B and the form (|12|) of the transformed potential 
into our condition (fTT]) : 

A (U )x - A u x u^ 1 U = -F A x u x u~ x - F A (u x u~ l ) x - A y u x u~ x - A (u x u" 1 )^ + 
+ {-A F x - A Uq + A u x u~ l F + A F A' 1 A x + A F u x u~ l + A y ) u x u~ l . 

This equation can be simplified by cancelling equal terms and factors on both sides, and 
by using ([9]). After regrouping terms, we arrive at the following form: 

(Uq) x - u x u~ x Uq + F x u x u~ x + F u xx u -1 + F u x (n" 1 )^ + u xy u^ 1 - 

—u x u^ 1 F u x u^ 1 + F u x u -1 u x n" 1 + u x (u~ 1 ) y + Uq u x u~ l = 0. (14) 
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We will now show that this equation can be integrated with respect to the variable x, 
which yields 

u- 1 (F u x + u y + U u) = C, (15) 

where C = C(y) is an arbitrary constant of integration. Let us now differentiate (|15p with 
respect to x and obtain equation (fT4"|) . Clearly, differentiation of the right hand side of 
(|15p gives zero, such that we only need to consider the right hand side: 

d x (vT 1 (F u x + u y + U u)j = (u" 1 )^ F u x + u~ x F x u x + u~ x F u xx + {vT v ) x u y + 

+ u~ l u xy + (u' l ) x Uq u + u~ 1 (U ) x u + u~ x Uq u x . 
Multiplication by u from the left and by u~ l from the right yields 

u d x _1 (F u x + u y + Uq uyj u~ l = u (u^ 1 ) x F u x u^ 1 + F x u x u~ l + F u xx u^ 1 + 
+ u (n -1 )^ u y u^ 1 + u xy u~ x + u (n _1 ) x U + (U ) x + Uq u x vT 1 . (16) 

Next, we make use of the following identity that holds for the derivative of an inverse 
matrix: 

(u~ 1 ) x = — u~ l u x u^ 1 . 

Clearly, this is also true for the partial derivative with respect to y, if x is replaced by y. 
Now we apply this rule to the derivative of the inverse matrices u~ l in ()16p and obtain 

u dxi^uT 1 (F u x + Uy + Uq uyj u^ 1 = —u x u" 1 F u x u~ x + F x u x u~ 1 + F u xx u^ 1 — 
- u x n" 1 Uy u^ 1 + u xy u" 1 - u x u~ x Uq + (Uq) x + Uq u x vT 1 . (17) 

Finally, on taking into account that 

- u x u^ 1 u y u^ 1 = u x (u~ 1 ) y , 

it is immediate to see that the left hand side of (|14p and the right hand side of (|17p 
coincide. This shows that our last condition (fTTJ can be integrated and turns into (fT5|) . 
which can be written in the form 

F u x + u y + U u - u C = 0. (18) 

Note that this equation turns into a matrix version of our Dirac equation Q if we choose 
C = 0. At this point the meaning of the matrix u becomes clear: each of its four columns 
contains a solution of the auxiliary equation (|18p . Thus, in order to perform our Dar- 
boux transformation ()22|) . we need a solution of the Dirac equation Q and four linearly 
independent solutions of the auxiliary equation (|18p . Now, we have completely resolved 
our conditions (|9j)- (|lip and we have found the intertwiner ([5]) that solves our intertwining 
relation ©: 

L = A (d x -u x u' 1 ) . (19) 

Consequently, for solutions ip and <f) of our Dirac equations ([3]) and (j3|), respectively, we 
have the desired mapping property Lip = <j). Next, we determine the transformed potential 
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U\ from its form (|12|) . On inserting the definition (|13|) . we obtain after simplification the 
following expression: 



U-i 



A (F x + Uq - u x u~ l F — F A' 1 A x + F u x u' 1 



A'' Ay) A 



-i 



(20) 



Now observe that the function A does not fulfill any constraint except that it must com- 
mute with F. We can therefore assume A to be a multiple of the 4x4 identity matrix 
I4, that is, A = a I4, where a is a constant that can serve e.g. for normalization purposes. 
In this case, the derivatives of A vanish and (|12|) simplifies to 



A (F x + Uq-u x u~ 1 F + F u x u' 1 ) A' 1 



U + F x + 



F, u x u 



-1 



(21) 



where [•, •] denotes the commutator. Hence, the solutions ip and <p of our Dirac equations 
(J3j) and (jlj), respectively, are interrelated by means of 



$ = Lip = A (ip x - u x u 1 ip) 
where the transformed potential U% is displayed by (|2ip . 



(22) 



3 The Dirac equation in cylindrical coordinates 

We will now review how the stationary Dirac equation in three dimensions can be written 
in cylindrical coordinates and how it can be reduced to a two-dimensional equation [8]. 
Afterwards, we show that our Darboux transformation is applicable. 

3.1 Introduction of cylindrical coordinates 

Let us start at the stationary Dirac equation in three dimensions: 

(i (ai d x + a 2 d y + a 3 d z ) - p m + E - Vb)) * = 0, (23) 

where the energy E is a real constant, Vq = Vo(x, y, z) stands for the potential, the positive 
constant m denotes the mass, and \P = VP(iE, y, z) represents the 4-component solution. 
Furthermore, recall that the matrices ay, j = 1,2,3, and (3 are given by 

where aj, j = 1, 2, 3, denote the usual Pauli matrices. Now we introduce the change from 
cartesian to cylindrical coordinates (x, y, z) 1— > (p, ip, z) by means of the usual relations 

x = p cos(ip) y = p sin(<£>) z = z. (25) 

Furthermore we suppose that the solution ^ of our Dirac equation (|23p depends on p and 
93, while its dependence on z can be separated off via 

^ = exp(z p z z) ip, (26) 
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where ip = ip(p, ip) and p z stands for the eigenvalue of the momentum operator p z = —i d z . 
Next, we substitute (|26p into our Dirac equation (|23p and rewrite the partial derivatives 
in cylindrical coordinates, which renders the Dirac equation in the following form [8]: 



Pi d p + ^ d v 



p z fc-P m + E-Vo \ V = 0, 



(27) 



where the matricial coefficients are given as follows: 
/ 



Pi 














\ exp(i <p) 



exp(— i if) 




exp(— i if) \ 
exp(z <f) 





(28) 



P2 



P3 



\ i exp(z if) 

( 1 







-i exp(— i ip) 






i exp(z ip) 





-i exp(— i ip) \ 






(29) 



\ 

-1 



1 

Vo 




-10 0/ 



(30) 



It is immediate to see that our Dirac equation in the form (|27p involves only the two 
variables p and ip, and that it is now a special case of our two-dimensional equation ([1]). 

3.2 The Darboux transformation in cylindrical coordinates 



Our goal is to make the Darboux transformation (|22|) applicable to the Dirac equation in 
cylindrical coordinates (f2T|) . such that we obtain solutions <f> and associated potentials V\ 
for the equation 



Pi d p + !j d v -p z P 3 - P m + E - Vt 



0. 



(31) 



The first step in applying our Darboux transformation is to convert the Dirac equation 
in cylindrical coordinates into the form ([3]). Since there are two terms containing partial 
derivatives on the left hand side of (|27p. we have two possibilities of arriving at the de- 
sired form of our Dirac equation. Consequently, we end up with two different Darboux 
transformations. 

First Darboux transformation. We multiply our equation (|27p by P^ 1 from the left: 

1 Pl l P2 



d p + 



d^-Vz Pi 1 Ps-Pi 1 P m + pY 1 (E-Vo)) V = 0, 



(32) 



Now we introduce the abbreviations 



Pi 1 P2 



U = - Pz P^ Ps ~ Pi 1 Pm + p^ 1 (E - Vo), 



(33) 
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which render equation (|32p in the form 



d p + 7 d v + U ] ip = 0. 



(34) 



On comparing this equation with our Dirac equation ([3]), we see that both coincide if in 
the latter equation we choose F = 7 and if we identify the variable x with ip and y with p, 
respectively. Thus, our Darboux transformation (|22p becomes applicable. Let u = u(p, <p) 
be a solution of the following auxiliary equation 



Up + 7 Utp + Uq u - u C 



0, 



(35) 



where C = C(p) is an arbitrary function. Clearly, this auxiliary equation is just (|18p with 
the same identifications of parameter and variables that has been done for (|34p. Now, the 
function <j) = cf)(p, <p) defined as 



4> = L tp 

is a solution of the Dirac equation 



- 1 



(36) 



9» + 7^ + f/i = 0, 



(37) 



where the transformed potential Ui is given by ()20p . which reads in the present case 



Ui = A(^ + U 



Up u 



-1 



7-7A l A^+~iu^u 1 - A 1 A p )A 1 



,-1 



ft" 1 /?2 ft" 1 /3 2 ^ _1 Ap ft" 1 /3 2 u v n- 1 



+ 



A~ x A p ) A- 1 



(38) 



where in the last step we made use of the fact that 7^ = 0. This expression can be 
simplified further, if we assume A to be a multiple of a constant a, that is, A = al^, where 
I4 stands for the 4x4 identity matrix. Recall that this assumption is valid, as the only 
constraint A must fulfill is to commute with 7. The transformed potential (|38p then takes 
the following form (recall 7^ = 0) 



Ui 



U + -fp + 



7, u v u 



U + 



1 

P 1 



ft 1 hi u <fi u 



-1 



(39) 



Thus, the solutions 1(1 and (j) of our Dirac equations (|34|) and (j37|) in cylindrical coordinates, 
respectively, are interrelated via our Darboux transformation (|36|) . where the auxiliary 
solution u and the transformed potential U\ fulfill ()35p and ()38p . respectively. Note that 
the potential difference (f39|) does not involve the actual potentials Vb and Vi, as they occur 
in the Dirac equation (|27p and in its Darboux-transformed counterpart (|3ip . The Dirac 
equation in its initial form (|27p obtained the form (|34p through a multiplication by P^ 1 ■ 
The inverse of this we now multiply by (|39p and make use of the relation (|33p : 



V: 



v + * 

P 



ft 1 ft,^ « 



(40) 



This is the final expression for the potential difference, where the potentials Vi and Vq are 
associated with the initial and the final Dirac equation (|27p and f)31 [) . respectively. 
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Second Darboux transformation. Let us now return to our Dirac equation it its form 
([27]) . and let us multiply this equation from the left by p P^ = p P^'- 



P P2 Pi d p + d v - Pz p fa Ps - P fa P m + p fa (E - V Q ) V = 0, 



(41) 



Similar to the previous case, we introduce the following abbreviations: 

7 = P fa Pi U = -p z p Pi fa- p fa P m + p p 2 (E - V ), 
which turn (|4ip into the equation 



(42) 



7 d p + &, + U ^ = 0. 



(43) 



Similar to our previous case (|34|) . we now compare our general Dirac equation ([3]) with 
(|43p . These equations coincide, if we have F = 7 and if we identify the variables p and (p 
with x and y, respectively. In this case we can apply our Darboux transformation 
the auxiliary function u = u(p, tp) for which must be a solution of the equation 



7 Up + u v + U u - u C 



0, 



(44) 



where C = C(ip) is an arbitrary function. Then, as in the previous case, the function 
(j) = 4>(p, (p) defined by 



cf) = L ip = A (Vp — Up u 1 i/A 
is a solution of the Dirac equation 



(45) 



7 dp + cL + Ui 



0. 



(46) 



The transformed potential U\ is given as the following particular case of (|20|) : 



-1 



A^Ap) A- 1 



Ui = A(^ p + U - Up u 1 7 - 7 A 1 A p + -f u p u 

= A (fa Px + Uq- p 2 up u' 1 P2P1-P P2 Pi A' 1 Ap + pfa fa up u^ 1 - A^Ap) , 

(47) 



In order to simplify this expression, let us suppose A = al^ for a constant a. The trans- 
formed potential (|47[) then takes the following form 



U + 7p + 



7, Up u 



U0 + P2 Pi+P 



P2 Pi, Up u 1 



(48) 



Hence, in the present case the solutions ij) and <p of our Dirac equations (|43p and (|46l) 
in cylindrical coordinates, respectively, are interrelated via our Darboux transformation 
(|45p . where the auxiliary function u and the transformed potential satisfy (|44p and (|47p. 
respectively. Let us finally remark that again relation (|48p can be converted into a relation 
between the actual potentials V$ and V\ in the Dirac equation (|27[) and its Darboux- 
transformed counterpart, respectively. To this end, note that the Dirac equation in its 
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initial form (|27p obtained the form (|43p through a multiplication by p The inverse of 
this we now multiply by (|48p and make use of the relation ()42p : 



= ^o + — + /3 2 [ft A,^^ 1 
P L 



(49) 



This is the final expression for the transformed potential V\ associated with the Dirac 
equation in its form (f3"Tj) . 



4 Examples 

We will now apply our Darboux transformation to the Dirac equation in cylindrical co- 
ordinates (|27p . for the sake of simplicity we consider the force- free case, that is, Vo = 0. 
To this end, let us first break up our Dirac equation into two components. We write its 
solution in the form tp = (ipi,ip2) T , where each of the functions tpi and tp2 stands for a 
two-component vector. Next, we need the following submatrices of the fij, j = 1,2,3, as 
given in (l28l)- (pOD : 



exp(— iip) \ _ / —i ex.p(—iip) \ __ / 1 
" 1 exp(iip) ) T2 ~ \ i exp(iyj) J T3 ~ \ -1 

(50) 

Now we are ready to decompose the Dirac equation (|27p into two equations, each of which 
has two components. 

n Wi) P H (^i)v> - Pz t 3 ipt + (E + m) tp 2 = 

P 

n {i>2) P H ($2)w - Pz t 3 ip 2 + (E - m) Vi = o. 

p 



Next, in order to apply our Darboux transformation (|22p . we need four linearly indepen- 
dent solutions of our auxiliary equation (|18p for Vq = 0, which will form the columns of 
our auxiliary matrix u. From now on we write this matrix u in 2 x 2 components as 

U\ u 2 

U 3 U 4 

For the sake of simplicity let us assume that the arbitrary matrix C = C(cp) that appears 
on the left hand side of our auxiliary equation ()18p . has the following form: 

c = ( V 2 c 4 °/ 2 ) • <W 

where each entry contains a 2 x 2 matrix, and both matrices C\ and C\ can depend on 
if. Let us now substitute (|5ip and the matrix u into the auxiliary equation (|18p . which in 
cylindrical coordinates then splits up into the following four 2x2 matrix components: 

^3 (52) 

p 

u\ (53) 

P 

■ U4 (54) 

P 

u 2 . (55) 

P 



(«a)p 


P 


(«3)<P 


~Pz 


r 3 u 3 + (E + m) u\ 


(ui) P 


+ ^ 
P 




-Pz 


r 3 m + (E - m) u 3 


(«4)p 


+ ^ 
P 




-Pz 


T 3 1i 4 + (E + m) u 2 


(«2)p 


+ ^ 
P 


(«a)»> 


-Pz 


t 3 u 2 + (E — m) u 4 
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We will now find two solutions of this system of equations, each of which leads to the 
generation of a new solvable Dirac potential in cylindrical coordinates. 



4.1 First example 

In our first example we take equations (|5'2p and (|53|) in the case E = m: 

n (u 3 ) p - p z t 3 u 3 + 2 m ui = ^ ^Ci u 3 - (ii 3 )^ (56) 

n {u\)p-p z r 3 u\ = ^ \ Ci ui - (ui^J. (57) 

In order to simplify these two equation further, let us now assume that C\ = (u 3 ) v = 0. 
We obtain the following equations: 

n (u 3 ) p - p z r 3 it 3 + 2 m u\ = 0, (58) 

7"2 

n (ui)p - p z t 3 m = (lilV (59) 

P 

Now we solve the first of these equations with respect to u\ and substitute the result into 
the second equation. The first equation yields 

Pz r 3 u 3 - ti (u 3 ) p 

ui = t, . (60) 

2 m 



the partial derivatives of which read 



Pz r 3 (u 3 ) p - ti (u 3 ) pp 



{m)p = ; ^ (6i) 

2 m 

= ■ (62 

2 m 

If we now replace ui in ()59p by inserting our findings (|60p . ()61|) and ()62p . we get 

- (u 3 ) P p -p 2 z u 3 = i {u 3 ) p (63) 

Let us suppose p z = 0, then this becomes 

- {u 3 )p P = i {u 3 ) p . 

A particular solution u 3 of this equation is given by 

n 3 = log(p) I 2 - 
Now that we have found u 3 , we can construct ui from relation (|60p : 

Tl 

u\ = ~= • 

2 m p 

It remains to find the matrix functions U2 and u^. To this end, we proceed in a similar 
way as done in the previous calculation. First we rewrite equations (|54p and ()55[) for the 

case E = —m, C4 = 0, (^2)^ = 0, and p z = 0: 

ti {u A ) p = [Uijtp, (64) 

p 

ti (n 2 ) p = 2 m Ui. (65) 
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The second equation can now be solved for u^. 

"4 = (u 2 )p- 

2 m 

The partial derivatives of this function read as follows: 

n 



(«4)p 



2 m 
2 m 



(«2) 
("2) 



PP 



P' 



Now we insert the function U4 and its partial derivatives into the remaining equation (|64p : 

u 2 = log(p) J 2 . 

As in the previous case, this is a particular solution. Now that we have found the four 
components of our matrix u, we can state this matrix in its explicit form: 

n 



( 



2 m p 
log(p) I 2 



log(p) I 2 



n 



(66) 



2 m p 



In order to find the transformed potential, the general form of which is displayed in (|47D . 
we need the partial derivative u p and the inverse u~ l of u: 



( 



n 



Imp 1 
P 



P 



2 m p 2 I 



n 



log 2 (p) + 



4m 2 p 2 

Multiplication of the latter two matrices gives 

/ /^ iog(p) L 



2 m p 
log(p) I2 



1 



log 2 (p) + 



1 



4 m 2 p 2 



n 



\ 2 



m /3 Z 



4 m 2 p 3 
1 + log(p) 



h 



log(p) I 2 
Imp 



n 

2 m p 2 
log(p) _ 



(67) 



1 + log(p) 
1 

4 m 2 p 3 



/ 2 



(68) 



We are now in position to construct the new solvable Dirac potential V\. To this end, 
we must evaluate relation (|4"9|) for the present case, which for the sake of convenience we 
repeat here, note that we have V\ = 0: 



- +/3 2 \fa Pi, Up u 1 



(69) 



We start by computing the first argument of the commutator in (j69j) , using the definitions 
and ([50]): 



Pi 



r 2 
r 2 



71 
71 



T 2 Ti 

t 2 n 



T 3 
T 3 



(70) 
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We will now compute the full commutator in (|69p . In order to avoid large expressions, we 
introduce the following abbreviations: 



A 



B 



1 



log 2 (p) + 



1 



4 m 2 p 2 



1 



log 2 (p) + 



log(l) _ 1 
p 4 m 2 p 3 



1 + log(p) 



2 m p 2 



4 m 2 p 2 

These settings render the product u p u~ l in the compact form 



A -B n 
B ti A 



Using this form of u„ u 1 we can finally evaluate the commutator, which reads 



[fo Pi, u p u 



-2 i 



r 3 
r 3 



B r 3 n 



A -B n 

B ti A 

-B T 3 T\ 





Next, we have to multiply this by (32 from the left: 



fa [fa Pi, Up u 



-li 



- 2i ( ° n 

-2 i B r 2 r 3 ri 





B T 3 Ti 

/ 2 
-h 



-B r 3 n 




The product r 2 r 3 t\ turns out to be 



(71) 



(72) 



(73) 

(74) 
(75) 



T~2 T 3 Ti 





exp (i ip) 



exp (— i 99) \ / 1 







exp [—% ip) 



— 1 J \ exp (i if) 



(76) 



Taking this into account, the commutator (|75p assumes the following final form, recall the 
definition of B as given in (|7ip : 



h [fo^upu- 1 ] = 2B 



(77) 



2 (1 + log(p)) 



2 m p 2 (log 2 (p) + T ^ I \ 

4 m (1 + log(p)) g 
4 m 2 p 2 log 2 (p) + 1 
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Now we take the previous result, combine it with 
/ 



, and substitute into 



1 



V\ 



exp 



H v) \ 



\ 



p 



o 



exp(i ip) 







exp(— i ip) 







P 



exp(z ip) 











+ 



/ 4 m (l+log(p)) 
' 4 m 2 p 2 log 2 (p)+l 



+ 



V 









/ 4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l 




- expfi <p) 
P 





4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l 








4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l 

— exp(— i if) 
P 








4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l 









1 



- exp(i y?) 

4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l 









4 m (l+log(p)) 
4 m 2 p 2 log 2 (p)+l / 

- exp(— i cp) ^ 



4 m (l+log(p)) 

4 m 2 p 2 log 2 (p)+l / 

(78) 



This is the final form of our Darboux-transformed potential V\. The corresponding solution 
(j) of the Dirac equation ()3ip can be constructed from the Darboux transformation (|45p. 

4.2 Second example 

We start again from the system of equations ([52]) - ([55]) . making the following assumptions: 

C 4 = (tfc)^ = = = 0. (79) 

Furthermore, we redefine E = e, together with the restriction < e < m. These settings 
render equations (|54|) and ()55p in the following form: 

(80) 
(81) 



n {ua) p + {e + m) U2 = 
ti (u2) P + (e - m) u 4 = 0. 

We can solve the first of these equations with respect to u<i'- 



U2 



e + m 



(82) 



On inserting this into equation (|8ip and introducing k 2 = e 2 — m 2 , we get 

{ua) pp = k 2 u A . (83) 
It is immediate to see that a particular pair of solutions to equations (|82p and (|83p 

T\ k cosh (A: p) 



u-2 



e + m 



U4 = smh(k p) I2 



(84) 
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It remains to find a pair of solutions for the equations (|52"j) and (follj) . To this end, ob- 
serve that under the settings ([79)) these equations become equivalent to (f54"j) and (f55j) . 
respectively, provided we set Ci=0. Thus, we will proceed as in the previous case, but 
introducing the setting E = —e. Note that this is allowed, since each pair of the functions 
(ui, U3) and (u2, ^4) must form a solution of our auxiliary equation (fTBj) . but not necessar- 
ily at the same energy E. Now, on employing (|79p . C\ = and = — e, we solve equation 
([52]) for U3 and substitute the result in (f53|) : 



«3 



n 



A 2 u\. 



e + m 

A particular solution of these two equations is given by 

Ti k cosh (A; p) 



(85) 



u\ = sinh(A p) I2 



e + m 



Now we have determined all four components of our matrix u, which reads in explicit 
form: 



V 



sinh(A; p) I2 
k Ti cosh(fc p) 



k Ti cosh(fe p) 

e + m 
sinh(A; p) I2 



£ + m 

We now perform the same steps as in the previous example. First we compute the inverse 
u^ 1 and the partial derivative u p of u, afterwards we use these results for evaluating the 
transformed potential. 



k cosh (k p) 
(e + m) 2 

/ 



+ sinh 2 (£i p) 



( 



\ 



sinh(A; p) I2 
k Ti cosh(fe p) 



\ 



cosh(fe p) I2 
k Ti sinh(fc p) 



£ + m 



£ + m 

k Ti sinh(fe p) 

£ + m 
cosh (A; p) I2 



k Ti cosh (A: p) 

£ + m 
sinh(A; p) I2 



The product of these two matrices appears in the transformed potential Vi, as given in 
(f4*9|) . and it reads 

_i 1 f m sinh(2 k p) I2 —k t\ cosh(2 k p) 

UpU ~ cosh 2 (A; p) + s ¥ k sinh 2 (A; p) \ k n cosh(2 hp) m sinh(2 k p) I 2 

Now, in order to evaluate (|49p . we introduce the following abbreviations 

m sinh(2 k p) 



A 



B 



k 



£ + m 



cosh 2 (A: p) 



£ + m 



sinh 2 (A p) I 2 



k cosh(2 k p) 



cosh 2 (A p) 



£ + m 



£ + m 

These settings render our matrix product u p u 



k 



sinh 2 (A p) 



-1 



u„ u 



-1 



A 

B Ti 



as given by ()68p . in the following form: 



-B Ti 
A 
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We now start to evaluate the commutator in (|49|) . Note that using our current abbrevia- 
tions A and B we obtain 

-B r 3 n 
B r 3 n 



[/3 2 /3i,«p u~ 



-2 j 



for the same reason as in (|73p . Furthermore, our calculations yield the following results: 



Pi [h Pi, Up u 



-2 i 



r 2 
r 2 



-2 i B T 2 T 3 Ti 





5 r 3 n 

h o 
o -J 2 



o 



(86) 
(87) 



where the matrix product t 2 t 3 ti simplifies to 

— exp (— i ip) 
exp (i if) 



T~2 T 3 Ti 





exp (i 



exp (— i ip) 




We substitute this result into (j87|) and get the commutator's final form: 



P2 [fo Pi, U p U- 1 ] 



2 B ft 



2 k cosh(2 k p) P 



k 



cosh 2 (A; p) H ^ — sinh 2 (/c p) 



e + m 
2 k 2 cosh(2 k p) P 



k 



—m + £ cosh(2 k p) 

On combining this with the definition (|28p and substituting into (|69p . we obtain 



/ 



Vi 



V 



p 







exp(i cp) 






exp(— i 99) 






exp(z 99) 





1 



/ 2 k 2 cosh(2 fc p) 
— m+e cosh(2 k p) 





+ 



V 








/ 2 k 2 cosh(2 fc p) 
—m+e cosh(2 fc p) 




- exp(i ip) 





2 fc 2 cosh(2 fc p) 
— m+e cosh(2 fc p) 








2 fc 2 cosh(2 fc p) 
—m+e cosh(2 fc p) 

- exp(— i (/?) 







exp(— i ip) 







2 fc 2 cosh(2 fc p) 
— m+e cosh(2 fc p) 




- exp(i cp) 
P 

2 fc 2 cosh(2 fc p) 
—m+e cosh(2 fc p) 





\ 



+ 








2 fc 2 cosh(2 fc p) 
—m+e cosh(2 fc p) / 

- exp(— i </?) ^ 





2 fc 2 cosh(2 fc p) 
—m+e cosh(2 fc p) y 
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This is the final form of our Darboux-transformed potential V\. As before, the corre- 
sponding solution (f> of the Dirac equation (|3ip can be constructed from the Darboux 
transformation (|45|) . 



4.3 Physical interpretation 

Let us now briefly discuss physical aspects of the transformed potentials that were obtained 
in the previous examples. To this end, recall the three-dimensional, stationary Dirac 
equation, minimally coupled to an electromagnetic field: 

ai (i d x + Ax) + a 2 (i d v + A 2 ) +a 3 (i d z + A 3 ) - j3 m + E - $ ^ = 0, 



where the electric charge was set to one and the definitions for ay, j = 1,2,3, and f3 
can be found in (|24p . Furthermore, <3? and (Ax, A 2 , A 3 ) denote the scalar and the vector 
component of the electromagnetic potential, respectively. We will now rewrite this Dirac 
equation in cylindrical coordinates, first observe the obvious fact that the equation can be 
written in the form 

f * ("i d x + a 2 d y + a 3 d^j + a t Ax + a 2 A 2 + a 3 A 3 - p m + E - ® \ 4» = 0. 



The differential operators translate into cylindrical coordinates as in (|27p . that is, we have 
after employing (|26p the following form of the Dirac equation: 

hid p + ^-d v> + a 1 A 1 +a 2 A 2 + a 3 A 3 -Pm + E-^)^ = Q, (90) 



where the definitions of Pi,f3 2 and f3 3 can be found in (|28|) - (|30|) . Before we see our two 
examples in detail, let us evaluate the term involving the vector potential components 
A 1: A 2 and A 3 : 



a\ A\ + a 2 A 2 + a 3 A 3 



( A 3 A x - i A 2 \ 

Ai + i A 2 -A 3 

A 3 Ai-i A 2 

\ Ax + i A 2 -A 3 J 



.(91) 



We will now see that the two transformed potentials constructed in the above examples can 
be written in the form (|9ip . such that the corresponding vector potential ($, Ax, A 2 , A 3 ) 
is determined. 



First example. We compare the matrix of the transformed potential Vx , as given in (JTHJ) , 
with the matrix (|9ip . Using Euler's formula, we first compare the off-diagonal elements, 
which leads to the following identification of the vector potential 

cos <A , sin 6 , , , 

$ = Ax = A 2 = A 3 = 0. (92 

P ' P 

It remains to accomodate the diagonal elements of our matrix (|78|) . This can be done by 
observing that the mass in the Dirac equation (|90p is a coefficient of the matrix j3 and can 
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thus be joined with the term (|77|) of the transformed potential V\. Thus, we introduce a 
position-dependent mass fa = fn{p) as follows: 



fa = m — 2 B 

= m 4m(l + log(p)) 
4 m 2 p 2 log 2 (p) + 1 

Hence, the transformed equation ([3"T]) with Vi given in ([75]) corresponds to a Dirac equa- 
tion with position-dependent mass (|93p . coupled minimally to a vector potential, the 
components of which are given in (|92j) . 



Second example. We proceed exactly as we did in the case of our first example. On 
comparing the matrices (|89p and (|9ip . we obtain the same identifications (|92|> as before. 
The position-dependent mass m = m(p) reads this time 

m = ra — 2 B 

2 k 2 cosh(2 fe p) 
— m + e cosh(2 fc p) 



Hence, the transformed equation ([31]) with V± given in ([89]) corresponds to a Dirac equa- 
tion with position-dependent mass (|94p . coupled minimally to a vector potential, the 
components of which are given in ()92|) . 

Note that in general case [18] is equivalent the following equation 

(M d x + N d y + E -V ) u = N uC(y), where E is arbitrary, (95) 

we can construct the exactly solvable transformed equation with any energy value E. 

( Ai \ ' 
A 2 
A 3 
\ A 4 / 

exactly solvable transformed equation (Vq — > V\) with matrix solutions v = {tp\,tp2i V3) V 9 4)) 
where (fi is solution to spinor transformed equation with E = Aj. 

Such as the initial equation with variable -E 1 , the solutions of transformed equation 
will with 'history'. We can consider appointed initial energy and appointed final energy. 
Also it is possible to construct solutions to transformed Dirac equation without 'history': 
v = Lv°, v° is matrix solution to initial Dirac equation. 

We hope that such results can be useful for study of electron transition, photon ab- 
sorption and creation. 



But if E is 4 x 4 diagonal matrix E = A 



we can construct the 



5 Conclusion 

The Darboux transformations for the general two-dimensional Dirac equation are con- 
structed. In fact, the Dirac equation in cylindrical coordinates is particular case of the 
general two-dimensional Dirac equation. Our calculations allow to construct the exactly 
solvable three-dimensional Dirac equation for a potential that depends on two variables. 
Exactly solvable potentials of the Dirac equation in cylindrical coordinates are generated 
by the Darboux transformation method, and the physical meaning of the transformed 
potentials is indicated. 
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